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1. Introduction 


The tsunami events that occurred over the last decades have caused an 
increase in public awareness and resulted in more research on the tsunami 
wave. Tsunami deposits play an important role not only in tsunami hazard 
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Tsunami inversion models attempt to link the basic information of the 
tsunami deposits with the overland flow characteristics^ There are three 
prominent inver sion m odels: M oore ’s advection model 


Souls by’s model 


Soulsbv et al. 


Moore et al. 


2007j and TsuSedMod model 


20071. 


Jaffe and Gelfenbuam 


2007] • It should be noted that all these models are based on different basic as¬ 


sumptions and employ different information from the deposits. For example, 
Moore’s advection model estimates tsunami flow magnitude by determining 


the combination of flow velocity and dept 
the sediment source to the deposition area 


l to move t 


le largest grain from 


Moore et al 


2007]. In this paper, 


we present a joint inversion framework (TSUFLIND), which combines these 
three models. TSUFLIND does not only couple all these three inversion 
models, b ut also conta ins a new method to calculate deposit characteristics 


Tang and Weis s 2014]. It also uses the calculated flow depth from Soulsby’s 


model to estimate a representative offshore tsunami wave amplitude. 


2. Theoretical Background 

2.1. Inversion Models Employed 

As mentioned above, there are three prominent tsunami deposition in¬ 
version models that will be used: Moore’s advection model, Soulsby’s model 
and TsuSedMod model. 


[Figure 1 about here.] 
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(a) Moore’s model:. 


Moore et al. 


2007] assumes that some grains in the sed¬ 


iment source do not move because the tsunami flow is not strong enough. 
Furthermore, it is assumed that most grains are transported in suspension. 
Based on these assumptions, the shear velocity is determined for the largest 
grain in the tsunami deposits. The law of the wall can be employed to find 
the shear stress, which is necessary to move the largest grain to get a flow 
velocity U. The following equation is used to determine deposition. 


h _ t _ l 
w, U 


( 1 ) 


in which w s is the settling velocity of the sediment grain, h is the water depth, 
l represents the horizontal distance a grain travels to be deposited. Because 
of the horizontal transport, this model is also referred as an advection model. 
This model was applied to d eposits forme d b y the 1929 Grand Banks 


tsunami, Newfoundland, Canada 


Moore et al 


2007]. In this application, it 


was estimated that the average flow depth was 2 .5 to 2.8 m, and the flow 


speed was 1.9 to 2.2 ms , which are the minima 


Moore et al. 


20071. 


(b) Soulsby’s model:. Soulsby’s model assumes that the water depth increases 
linearly between 0 and pT and decreases from pT to T for any given locations. 
T is the inundation time and p is a constant related to run-up time, which 
is between 0 and 1. H = H 0 + Ah is the maximum flow depth at a given 
location during tsunami inundation and decreases toward the inundation 
limit, Ho denotes the maximum water depth at the shoreline, Ah denotes 
the depth increment due to tsunami: 
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Ah = 


1{R Z ~ H 0 

mR? 


l 

m 


( 2 ) 


where m is the slope and R z represents the vertical inundation limit. The 
thickness of the deposit for grain size i at the shoreline: 


Ai) _ 
SO — 


C-Q Ws Id 


-(1 + 7) 


( 3 ) 


(1 -p)p s (l + ccW) 

where a ^ = Ws H ^ d , denotes the settling velocity for grain size i, T d = (1 — 
y)T is the deposition time. Cq is the depth averaged sediment concentration 
for grain size i and p is the porosity. The sediment thickness for grain size i 
linearly decreases with distance from the shoreline: 


. CS°(1--^) x<Rf 
C w (x) = 1 

' 0 x>R® 


( 4 ) 


where r!'s } is the distance between sediment extend and the shoreline for 
grain size i 


Soulsby et al 


2007] 


(c) TsuSedMod:. 


Jaffe and Gelfenbuam 


2007| developed an inversion model 


based on sediment deposited from suspension. There are several basic as¬ 
sumptions in TsuSedMod: (1) sediment is transported in suspension and 
deposited when steady and uniform tsunami flow slows down; (2) suspended 
sediment concentration is distributed in an equilibrium profile; (3) there is 
no erosion caused by return flow. The model iteratively adjusts the sediment 


source and the shear velocity to match the sediment 


grain-size c 


Jaffe et al. 


and thickness of suspension-grading sediment layers 
For the grain size i , the sediment thickness ApW is given by: 


istri 


mtions 


2011 


2012 ], 
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i r H ( x ) 

Ar/ (i) = -- / C®( z )dz (5) 

(1 — P) Jo 

where C^\z) is the sediment concentration profile of grain size i. After 
determining the shear velocity, the flow speed profile is calculated by : 


U{z) = 


rZ ul 




dz 


( 6 ) 


where z 0 is the bottom rough ness from 


eddy viscosity profile from 
The TsuSedMod model 


MacWilliams 


Gelfenbaum and Smith 


1986] 


2004 1 and K(z) is the 


l as been app 


Jaffe and Gelfenbuam 

2007, 

Spiske e 

al. 

2010; 

Jaffe et al. 

2011, 

2012 

and 

two paleotsunami cases 

Witter et al. 

2012 

; Soiske et ah 

201 

3 . For the 2009 


ied to four modern tsunami cases 


tsunami near Satitoa, Samoa, the flow speed estimated from TsuSedMod at 
three locations (100, 170 and 240 meters inland) were 3.6 to 3.8 ms" 1 (bottom 
layer/earlier wave) and 4.1 to 4.4 ms" 1 (top layer/later wave). These results 
are consistent with the 3 to 8 ms" 1 flow speed from the boulder transport 


inverse model 


Jaffe et al, 


2011 1 


2.2. Sedimentation model 

The method used to calculate the sediment concentration of the sediment 
source in TSUFLIND is similar to the one presented in 


Madsen et al. 


1993], 


The grain-size distribution of the sediment source is characterized by D 50 , 
the largest grain and the smallest grain. 

When the entire tsunami deposit at a given location is considered, re¬ 
suspension sediment flux can be neglected and Soulsby's model is applied. 
However, if the individual layer in the tsunami deposit is considered, intense 
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turbulent mixing cannot be ignored. Therefore resuspension has to be taken 
into account. The generation of each individual portion of the tsunami sedi¬ 
ment based on flow condition is the fundamental part to reconstruct tsunami 
deposits. For the entire deposit, the basic process is to calculate sediment 
thickness Q l \ x ) f° r each grain size at each point along the slope by using 
Eq ([3j) and (J4J) from Soulsby’s model. We assume that the depth averaged 
sediment concentration Co in Eq ([3j) is the reference sediment concentration 
C r here. The refere nce co ncentration is calculated for a given flow condition 


with 


Madsen et al. 


1993J: 


r< „ = ft(l-p)/ (i) 5 w m 

1 + /?„S«> [l 

where (3o is the resuspension coefficient, /W is a fraction of the sediment of 

grain size i. S W is the normalized excess shear stress given by 


= 


Tb-T c 

~U) 


(0 


T b > T^r 


0 T b < Ter 

(0 


( 8 ) 


where r b is the bed shear stress and Ter is the critical shear stress of the 


initial sediment movement for grain size i 


Madsen et al. 


1993 1. 


For a given location x, the grain-size distribution for the entire tsunami 
deposit is given by: 


/ 


(0 


C (i) ( x ) ■ 

Ef =0 C«(*)’* 


1,2,3, 


,N 


(9) 


where is the percentage of grain size i in the entire sediment, (^(x) is 
sediment thickness of grain size i and E total deposit thickness 

for all grain sizes. N is the number of grain size classes. 
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The tsunami deposit characteristics are reconstructed by matching sedi¬ 
ment thickness and grain-size distribution with held data. In order to recon¬ 
struct deposit details, the sediment concentration cannot be depth averaged 
and is described as a Rouse-type suspended sediment concentration prof ile. 


In this framework, we use the method from 


Jaffe and Gelfenbuam 


3 


2007 to 


calculate the suspended sediment concentration profile. It is efficient to re¬ 
construct the deposit by calculating times of deposition. The deposition time 
of suspended sediment is calculated by: 


j-0) _ 


w, 


(0 


( 10 ) 


in which i!'- } is the deposition time for grain size i sediment at elevation Zj. 
The amount of sediment settling in each grain size class for a given elevation 
is tracked as sediment thickness increment: 


cf = ’ 


1 ~P 


( 11 ) 


Cf' 1 is from the suspended sediment profile 


Jaffe and Gelfenbuam 


2007J, C 


(0 


is the thickness of the sediment increment of the same grain size i at elevation 
Zj and deposited at time tS l \ The deposition time and corresponding sedi¬ 
ment thickness increment are ordered from shortest to longest. If there are 
multiple layers in the tsunami sediment, we can compute the grain-size dis¬ 
tribution for each layer separately based on the depositional temporal order 
of the sediment thickness increments by: 


Ai) _ 
Jk ~ 


y M c 

2^.7=0 v 


(0 




Z-^i=0 \ Z-/7=0 3 ? 


M Ai) 


j i — 1)2,3,..., N\ j — 1, 2, 3, 


M 


( 12 ) 
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where is the sediment fraction of grain size i in layer k. o is total 
sediment thickness with the same grain size i in sediment layer k. Index 
j is used to mark the original location of sediment in the water column. 
YliLo ( S/=o ) i s the total thickness of this sediment layer which contains 


all grain size classes. In TSUFLIND, the calcula tion of tsunami flow condition 


will use the same method as TsuSedMod model 


Jaffe and Gelfenbuam 


2007|. 


2.3. Result Evaluation 

We employ the second norm to quantify the error between model and 
observed results as a control of the iterative procedure. The second norm of 
error for layer k is given by: 

, ,/e,". (f$-fP) 2 

Li ~ V- N - 

fm and fo ' 1 are the modeled and observed percentages for each grain size 
class i. With the help of Lk , we compute the average second norm value for 
a location with: 



£ = < 14 > 

k =1 

We define L < 5% as a good simulation. For the tsunami sediment thickness 
simulation, we employ the same process. The second norm value of error for 
thickness between the model result and the held observation is given by: 


Lth — 




2 


(15) 











where th m and thf are the modeled and observed thicknesses for each sample 
location, Q is the number of sample locations. As there is only a limited 
number of tsunami deposit samples for the test case applied here, we use 
10% as the threshold value. 

2-4- Offshore Wave Characteristics and Flooding 

In order to estimate a representative offshore tsunami amplitude, we re¬ 
late the water volume calculated from Sousby’s model at maximum inun¬ 
dation with the volume calculated by numerically solving the shallow water 
equation. We carry out a parameter study by varying the slope (m) and the 
offshore wave amplitude (£). For more details about the parameter study 
and employed numerical model, we refer to Appendix II. The water depth 
computed from Soulsby’s model is used to calculate the volume of the in¬ 
undation water. With the help of numerical simulations (Appendix II), we 
derived the following formulation: 

£ =Ai T A2 ■ V T A3 • m T A4 • V T A5 • m • V T Ag • 

(16) 

T A7 • V '' T As ■ V“ • m T Ag • V • mf T A^q • mf 

Where £ is offshore wave height, V is the water volume that covers the land 
at maximum inundation, m is the slope of beach profile. These constants 
A in Eq 16 are Ai = 5.06, A2 = 2.93, A3 = —0.28, A4 = 0.51, A5 = —3.04, 
A 6 = 0.0014, A 7 = 0.027, A 8 = -0.011, A 9 = 0.051, Ai 0 = 0.053. 

[Figure 2 about here.] 
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2.5. Inversion Framework and Coupling 

We use the information from all three models as different components in 
this joint inversion framework. The steady flow condition that is presented in 
all models, is also presented in TSUFLIND and represents the most simplify¬ 
ing assumption. The inputs to TSUFLIND are the sediment characteristics 
for different sampling locations along a slope. However, it should be noted 
that the inversion of the flow conditions is carried out for each sample loca¬ 
tion individually. TSUFLIND uses components from Moore model, Soulsby’s 
model and TsuSedMod model to adjust the sediment source grain-size dis¬ 
tribution, the sediment source concentration and the average flow velocity 
to simulate tsunami sediment thickness and grain-size distribution along the 
slope in the deposition zone. If needed, the representative offshore wave am¬ 
plitude can be computed. Figure 2 depicts the flowchart outlining how the 
joint inversion model works. 

The information needed for a successful inversion includes the grain-size 
distribution, sediment thickness as well as the information of the slope along 
which the tsunami sediments were sampled. In the inversion framework, 
Moore’s advection model is employed to calculate the initial flow speed. Be¬ 
cause the Moore’s model uses the actual data from the measured sediment 
distribution, it reduces the number of iterations significantly. The reservo ir 


of sediments in the water column is calculated by following 


Madsen et al. 


1993] . and it is assumed that all grain-size distributions can be described 


with log-normal distributions. The iteration begins with computing the in¬ 
undation with the help of Soulsby’s model, and the initial estimates of the 
flow conditions are from the Moore’s advection model. The results of this 
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step are the local flow depth and the entire sediment thickness at each sam¬ 
ple location. Our sediment formation module calculates the characteristics of 
the deposited sediments. The iterations are controlled by the norm of error 
between the simulated and observed deposits and stop after the predefined 
threshold is met. As the model outputs, we can estimate the flow speed, 
depth and Froude number along the slope. If needed, a range of offshore 
reference wave amplitudes can also be computed. 


3. Application and Example 

|Figure 3 about here. 

3.1. Field Observation and Data 
We employ the field data 


Ba hlburg and Weissll2007| from the 2004 Indian 


Ocean tsunami to demonstrate the capabilities of our framework (Figure 3). 


These samples come from th e coas 


puram, India. Bahlbur g and Weiss 


al are a in the vicinity of Ranganatha- 


2007] identify sediment layers formed 


by the tsunami in this cross section and described grain-size distributions for 
each layer. There are some grass runners on the top of the tsunami sediment, 
which indicate the return flow direction and the erosion caused by the return 
flow. Most grain-size distributions of the sediment layers in the test case 
are unimodal (Figure 3b). Tsunami deposits in this cross section are usually 
well sorted, and the mean grain size is between 0.5 and 1.5 in 6 scale, which 


corres ponds to medium and coarse sand. Furthermore, Bahlbu rg an d Weiss 


2007] observe that the mean grain size is upward and landward fining in this 


cross section. Several sedimentary data, such as the deposit thickness and 
the grain-size distribution, will be used as input to TSUFLIND. The range of 
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the grain sizes in the sediment source is based on all held samples collected 
during the held survey. Flow depth in this model will take full use of both 
the held observations and the model results from Soulsby’s model. 

3.2. Sedimentary Simulation Results 

TSUFLIND hrst simulates tsunami deposit thickness (Figure 3c). In the 
test case, the largest observed thickness is about 0.22 meters at 120 meters 
inland. For the hrst 100 meters in this cross section, the simulated thickness 
from TSUFLIND is obviously larger than the held measurement. After 200 
meters inland, the simulated thicknesses decrease quickly and generally ht 
with the held measurement. 

TSUFLIND reconstructs sediment grain-size distributions for both the 
entire tsunami deposit and several vertical intervals at any given sample 
locations. The error of the entire tsunami sediment grain-size distribution 
in this test case is from 0.38% to 1.54%, which can be considered good 
simulation results. The error is less than 1.0% from 120 meters to 160 meters 
inland and then increases to 1.5% after 160 meters inland. We use four 
sediment samples to calculate grain-size distributions (Figure 3d I — IV , 
response to 120m, 160m, 177m and 207m from shoreline). Beyond 160 meters 
inland, there are fewer coarse grains and more hne grains in the simulated 
grain-size distribution than the held measurement (Figure 3d, /, III and 
IV). 

In order to study how the grain-size distribution changes in the vertical 
direction, we employ the new sediment formation module to simulate tsunami 
deposit grading. Figure 3a shows grain-size distribution for several vertical 
intervals at four different study locations. The grain size for these recon- 
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struction results ranges from 0 to 6 in 0 scale, so the portions of simulation 
results that are outside of this window are not plotted in Figure 3a. The 
number of vertical intervals separated decrease away from the shoreline as 
the sediment thickness decreases to zero. These simulation results show some 
features similar to the field observation, such as fining inland and upward. 

Based on the grain-size distribution for each vertical interval shown in 
Figure 3a, the sediment parameters used to describe the deposits can be 
calculated. Taking the sample from 120 meters inland as an example, we 
calculate the mean grain size, kurtosis, skewness and sorting factor for each 
interval. In the bottom several centimeters, the mean grain size does not 
change significantly and is about 1.2 in 0 scale. Then the mean grain size 
decreases upward to 2.2 in 0 scale. The change in kurtosis is about 0.8 to 1.1 
in this sample, which indicates the grain-size distribution has a wider grain 
size range than a normal distribution. Sediment simulation results in this 
example also show that tsunami sediment changes from moderate sorting at 
the bottom to well sorting at the top where fine grains become dominant. 
Also the grain-size distribution is positively skewed which indicates that the 
distribution is skewed to fine grains. 

3.3. Hydrodynamic Inversion Results 

[Figure 4 about here.] 

After reconstructing grain-size distributions, TSUFLIND calculates the 
flow speed and Froude number at the sample locations. In the test case, Fig¬ 
ure 4a and 4b show the flow speed and Froude number distribution along the 
slope. The average flow speed decreases from 4.7 ms -1 at 150 meters inland 
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to 3.3 ms -1 at 350 meters inland. The Froude number, which is around 0.9, 
does not change significantly in space. Also the range of possible velocities 
and Froude numbers decreases from 150 meters to 350 meters inland. 


The flow speed profile shape is influenced by the eddy viscosity profile 
and shear velocity. The eddy viscosity profile parameterizes the vertical 


distribution of turbulent stress. TSUF 
profile based on laboratory data from 


IND foll ows t he flow eddy viscosity 


Gelfenbaum and Smith 


1986 j. The 


flow speed gets the largest value on the water surface and decreases toward 
the sediment bed. TSUFLIND only gives the depth-averaged values as final 
results. 

TSUFLIND computes the water surface profile to estimate the water 
volume when the tsunami wave reaches the maximum inundation. With the 
help of Eq (16), the wave amplitude can be estimated based on the slope (m) 
and the water volume (D). For the Eastern India case, the wave amplitudes 
range from 5 to 7 meters and the wavelength is close to 50 km. 


4. Discussion 

4-1. Interpretation of test case results 

With the help of the presented model, we can reproduce tsunami sedi¬ 
ments as well as infer the flow condition that generated them based on ob¬ 
servations and analytical results of existing tsunami deposits. Figure 3 sum¬ 
marizes the results of our simulation for the tsunami deposits. The apparent 
difference of the deposit thickness between model results and observations for 
distances smaller than 120 m from the shoreline can be explained by strong 
return flow or large velocities from subsequent waves with small inundation. 
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For distances from the coastline larger than 120 meters, the deposit thickness 
simulation results matched well with the observations. However, the obser¬ 
vations are slightly larger due to the presence of topography change that may 
slow down the flow (Fig 3c and d). The finer grain-sizes contain the largest 
error between observation and model results is likely that the topography 
change is the main source of the error. However, the difference could also be 
a part of the model uncertainty. 

The calculated mean speed decreases from 4.7 ms -1 to 3.3 ms -1 along 
the studied section. The speed decreases continuously shown as Fig 4a, the 
Froude number increases and then decreases (Figure 4b). The mean Froude 
number is around 0.9 for this test case. As the change rate of water depth is 
a constant based on Sousby’s model, the mean Froude number changes only 
depends on the velocity changes. At first, the flow speed decreases slower 
than the water depth, so the mean Froude number increases in this area (150 
m to 300m). After some point around 300 meters, the flow decelerates quickly 
and causes the Froude number to decrease. The flow speed and Froude num¬ 
ber results from TSUFLIND are shown as ranges of possible values with 
uncertainties (Figure 4a and 4b). In Figure 4, the gray area with solid line 
boundaries is the estimated flow speed or Froude number with maximum 
and minimum possible values. The dashed line is the mean value of the esti¬ 
mated flow speed or Froude number. The ranges of the speed and the Froude 
number decrease from 150m to 350m, which indicates the uncertainties de¬ 
crease towards the sample location close to the landward sediment pinch-out. 
It is possible that the tsunami deposits near the maximum run-up position 
become thin, well-sorted and fine-grained containing less information about 
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the flow condition. Tsunami wave amplitudes calculated by TSUFLIND are 
usually larger than real amplitudes, because the mathematical relationship 
(Eq 16) is from frictionless shallow water equations that ignore the effect of 
friction. 

4-2. Model limitation and improvement 

In this study, we combine three tsunami inversion models to simulate 
tsunami deposit and estimate tsunami flow parameters. All three of these 
models are based on some model-specified basic assumptions. A significant 
assumption of TSUFLIND is that the sediment transport and the deposition 
process during a tsunami are considered uniform in space and time. Conse¬ 
quently, the deposit conies from both horizontal convergence and suspension 
settling. TSUFLIND combines Sousby’s model and TsuSedMod to simulate 
these two processes. This combination greatly improves the grain-size dis¬ 
tribution simulation results. However, when the tsunami flow decelerates 
rapidly because of bathymetry changes or any other reasons, some part of 
the sediment deposited before will be eroded again. If the flow is strong 
enough, a significant part of tsunami deposit may be eroded, just like the 
result shown in Figure 3c from shoreline to 100 meters in land. As a result, 
the tsunami speed calculated by TSUFLIND will be underestimated. 

Another significant assumption of TSUFLIND is that most of tsunami 
deposits is transported by the suspension load and ignores the contribu¬ 
tion of the bed load. This assumption will overestimate the tsunami flow 
speed and increase the percentage of coarse grains in the grain-size distribu¬ 
tion. TSUFLIND is not applicable where bed load is the dominant sediment 
transport type. In most tsunami cases, tsunami deposits include both bed 
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load and suspended load. In order to reduce the effect of bed load, only 
the suspension-grading part can be considered to estimate the flow speed. 
Suspension-grading in sediment is difficult, to identify only based on the grain- 
size distribution. The simulated suspended-grading part has some difference 
with held observations, which also introduces uncertainties to the how speed 
calculation. 

TSUFLIND has two aspects that require improvement: the applicable 
how condition and the accuracy of hnal outputs including sediment simu¬ 
lation and how speed calculation. The improvement of the applicable area 
can be made by employing other tsunami propagation models instead of 
Sousby’s model, which can deal with non-uniform and unsteady hows. At 
the same time, the new model needs to consider both the suspension load 
and bed load. Also a new method for combining forward and inverse mod¬ 
eling will hold great promise for deciphering quantitative information from 
tsunami deposits and decreasing the u ncerta inties in tsunam i sediment trans¬ 


port simulation and inversion results 


Sug awar a et al. 


20141. 


5. Conclusion 

Modeling the tsunami sediment deposition processes and estimating tsunami 
how parameters will greatly improve not only the understanding of deposition 
from tsunami but also the risk assessment for extreme high-energy events. 
We presented a joint inversion model for tsunami deposit simulation and how 
condition estimations in this paper. 

The combination of different inversion models allows the computation of a 
wide range of tsunami wave impacts or characteristics, ranging from sediment 


17 





thickness, grain size distribution to flow speed and wave amplitude. However, 
these results are restricted by the flow condition. For instance the flow speed 
of non-uniform or unsteady flow cannot be inversed yet. Also if there were 
strong return flows or subsequent waves, the results from TSUFLIND will 
contain some uncertainties. From a general point of view, with a simple 
bathymetry, the tsunami and paleotsunami as well as other extreme events 
can be understood with the help of this joint inversion framework. More 
research needs to be done to improve this joint inversion framework presented 
to quantify and reduce the uncertainties in the inversion results and expand 
applicable conditions. 

6. Appendix I 


Table 1: Symbols List 


Symbol 

Unit 

Description 

h 

m 

Water depth 

W s 

ms -1 

Settling velocity of the sediment grain 

l 

111 

Horizontal distance a grain travels to be deposited 

7 

1 

Run-up time constant 

T 

s 

Inundation time 

H 

111 

Maximum flow depth at a given location during tsunami inundation 

Ah 

m 

Water depth increment due to tsunami 

m 

i 

Slope 

Rz 

m 

Vertical water inundation limit 

H 0 

m 

Maximum water depth at the shoreline 

Ai) 

to 

m 

Thickness of the deposit for grain size i at the shoreline 

T d 

s 

Deposition time 

°0 

m 3 /m 3 

Depth averaged sediment concentration for grain size i 
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Table [T] Continued: Symbols List 


p 

1 

Porosity 

£(*) 

m 

Sediment thickness for grain size i 

R ( s l) 

m 

Distance between sediment extend and the shoreline for grain size i 

A r)W 

m 

Sediment thickness increment for grain size i 

C7W(z) 

m 3 /m 3 

Sediment concentration profile for grain size i 


m 

Bottom roughness 

*(*) 

kgm _1 s _1 

Eddy viscosity profile 

Cr 

m 3 /m 3 

Reference sediment concentration 

A> 

1 

Resuspension coefficient 

/« 

% 

Percentage of the sediment of grain size i 

S« 

1 

Normalized excess shear stress 

n 

kgm -1 s -2 

Bed shear stress 

T « 

1 cr 

kgm -1 s -2 

Critical shear stress of the initial sediment movement for grain size i 

N 

1 

Number of grain size classes 

z j 

m 

Sediment original elevation 

A*) 

l 3 

s 

Deposition time for grain size i sediment at elevation Zj 

Lk 

1 

Second norm of error for layer k’s grain-size distribution 

/■(*) 

Jm 

% 

Modeled percentages for grain size class i 

Ai) 

Jo 

% 

Observed percentages for grain size class i 

L 

1 

Average second norm value of grain-size distribution for a location 

Lth 

1 

Second Norm value of thickness between the model result and the field observation 

thm 

m 

Modeled thicknesses 

thf 

m 

Observed thicknesses 

Q 

1 

Number of sample locations 

£ 

m 

Offshore wave amplitude 

V 

m 3 

Water volume that covers the land at maximum inundation 

Rw 

m 

Distance of water run-up to shoreline 

v(x,t) 

m 

Solitary wave form 

C 

ms' 1 

Wave celerity 
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Table [T] Continued: Symbols List 


u(x, t) 

ms 1 

Horizontal velocity in shallow water equations 

9 

_2 

ms 

Gravitational constant 

B(x) 

m 

Bed topography function 

d 

m 

Water depth of continental shelf 


7. Appendix II 

TSUFLIND implements a simplified method to estimate the representa¬ 
tive offshore tsunami wave amplitude. First of all, the water volume on the 
land due to tsunami wave is calculated by: 

r Rw 

V — / h(x)dx (17) 

Jo 

where V is the water volume, R w is the distance of run-up to shoreline, h(x) 
is water depth distribution function on land. To simplify this problem, we 
assume tsunami wave is a solitary wave. The solitary wave form is given as 
a function of distance x and time t by 

rj(x,t) = £sech 2 (k(x — Ct )) (18) 

where 

k =M < i9 > 

£ is the wave amplitude and d is the water depth of continental shelf, which 
is assumed as 500 meters. C is the wave celerity which is expressed as: 


C=y/g(h+H) 


( 20 ) 
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Initial velocity in shallow water equation code is set as: 


vo (x,t) = \/g/h O i((,,fco) 


( 21 ) 


We calculate the water volume when tsunami wave got the maximum run¬ 
up based on water distribution function h(x). Water depth function h(x) 
comes from a shallow water equations code. The shallow water equations 
code used here originally is designed for studying the propagation and runup 
of solitary wave by using a high-resolution finite volume method to solve 


following equations 


Delis et al. 


2008|: 


dh d{uh ) 


dt 


dx 


= 0 


( 22 ) 


diuh) d ,, n 1 , dB 

m_i + _ (te2 + _ 3h 2 ) = _ 9ft _ 


(23) 


where h(x,t) is the flow depth, u(x,t) is the horizontal velocity, g is the 
gravitational constant, B(x) is the bed topography function. 

In this code, a conservative form of the nonlinear shallow water equations 
with source term is solved by using a high-resolution Godun ov-typ e explicit 
scheme with Roe’s approximate Riemamn solver 


Delis et al 


20 08 ]. In order 


to get the mathematical relationship between the maximum water volume 
(V), slope (m) and initial wave amplitude (£), we design a parameter study 
by varying slope and wave amplitude to calculate the water volume. Finally, 
we use curve fitting methods to get the mathematical relationship based on 
parameter study data set. 
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Offshore 


Erosion zone 



Figure 1: Conceptual model of TSUFLIND with definition of the terminology 
used later in the paper. For more symbols used in this paper see Appendix I. 
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Figure 2: Flowchart for TSUFLIND’s iterative scheme to simulate tsunami de¬ 
posit and estimate tsunami flow condition. 
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Figure 3: TSUFLIND simulation results and field measurement at Ranganatha- 
puram, India, a: Vertical grading in grain size distribution (blue line) and mean 
grain size (red line) for four sampled locations (120m, 160m, 177m and 207m); b: 
the entire tsunami deposit grain-size distributions used as inputs to TSUFLIND 
(red points) and model result outputs from TSUFLIND (green line); c: tsunami 
deposit thickness field measurements (red points) and simulation results from TSU¬ 
FLIND (green line); d: topography, wave run up and sample locations for test case 
(I : 120m; II : 160m; III: 177m; IV : 207m). 
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Figure 4: The estimated flow speeds and Froude numbers from TSUFLIND. a: 
Tsunami flow speed estimates are indicated by the gray area with the boundaries 
of maximum and minimum possible speeds. The dashed line is the average value 
of estimated flow speeds, b: Froude number estimates are indicated by the gray 
area in this figure with the maximum and minimum possible values. The dashed 
line is the average value of Froude number. 
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